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Abstract 

An iterative method for recovering the bulk information in asymptoticahy 
AdS spacetimes is presented. We consider zero energy spacehke geodesies and 
their relation to the entanglement entropy in three dimensions to determine the 
metric in certain symmetric cases. A number of comparisons are made with 
an alternative extraction method presented in arXiv:hep-th/0609202, and the 
two methods are then combined to allow metric recovery in the most general 
type of static, spherically symmetric setups. We conclude by extracting the 
mass and density profiles for a toy model example of a gas of radiation in 
(2+l)-dimensional AdS. 

1 Introduction 

One of the original applications of the holographic principle was in relating the 
entropy of a black hole to the area of its horizon [U [2]; since then a variety of 
authors have continued to explore relationships between bulk and boundary physics 
via holography, most notably via the AdS/CFT correspondence of Maldacena [3]. 
The idea of entropy being linked with an area rather than a volume (as one naturally 
expects from thermodynamics) is not, however, restricted to the case of black holes. 

Recently, a proposal was put forward by Ryu and Takayanagi [H [5] relating the 
entanglement entropy of a subsystem in a CFT to the area of a minimal surface 
in the bulk. This has been investigated further in a number of subsequent papers, 
such as El El E] where a number of related issues are explored. One avenue of 
interest leading from this proposal is the question of whether we can take this link 
between the entanglement entropy and minimal surface area, and devise a method 
to efficiently extract the bulk physics from the field theory information. 

In (2+1) dimensions, the area of the minimal surface in question corresponds to 
the length of a static spacelike geodesic connecting the two endpoints of the region 
A through the bulk, as illustrated in figure [2j It is this observation that leads to 
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comparisons with a method of extracting the bulk metric given in [lOj, where the 
relation between singularities in correlation functions in the CFT and null geodesies 
(see [TT] for details) was used to iteratively recover the bulk metric in certain asymp- 
totically AdS spacetimes. In this paper we devise a similar method for extracting 
the bulk metric, using instead the relationship of Ryu and Takayanagi between the 
entanglement entropy and the length of the relevant spacelike geodesic. Interest- 
ingly, we find that after plotting the proper length against the angular separation 
of the endpoints, see figure El the gradient dC/dcf) immediately yields the angu- 
lar momentum of the corresponding static spacelike geodesic. This simple relation 
then allows the minimum radius of the geodesic to be determined, and by working 
iteratively from large r, one can reconstruct the metric function of the bulk. 

After describing the method and giving some examples of its application in 
practice, we then make a number of comparisons between this and the method 
of [TO] (which is briefly reviewed in section 14. ip . Most crucially, the two methods 
involve different ways of probing the bulk (as they involve different types of geodesic 
path), and whilst they appear computationally quite similar, this difference allows 
the method presented here to probe more fully a greater range of asymptotically AdS 
spacetimes. This is a consequence of the fact that in singular spacetimes, and those 
with a significant deviation from pure AdS, the effective potential for the null paths 
can become non-monotonic, resulting in geodesies which go into unstable orbits, see 
figure [T2I This local maximum in the potential results in a finite range of radii 
which cannot be effectively probed by the null geodesies, and information about the 
bulk cannot be extracted; one does not encounter this problem when probing with 
static spacelike geodesies, provided the metric function is non-singular. Despite this 
advantage, one cannot use either method individually to extract information from 
the most general static, spherically symmetric spacetimes (those with a metric of 
the form of ()26p ). as neither can provide enough data with which to fully determine 
the metric; the null geodesies are not sensitive to the overall conformal factor of 
the metric, and the static spacelike geodesies cannot probe the timelike part. One 
can, however, use them in conjunction in order to do so. We thus conclude by 
proposing a combination of the two approaches such that the bulk information can 
be recovered, and give firstly an example demonstrating the ease with which it can 
be done, followed by a toy model setup of a gas of radiation (a "star") in AdS^^. 
We demonstrate how it is possible to determine both the star's mass and density 
profiles from our estimates of the metric functions. 

The outline of the paper is as follows: Section [2] contains background material on 
asymptotically AdS spacetimes and geodesic paths, and introduces the entanglement 
entropy relation from [3j . Section [3] develops the method for iteratively extracting 
the bulk metric, the full details of which are given in Appendix A, comments on 
the validity of the solutions, and goes on to give examples. In Section HI after a 
review of the null geodesic approach from [IQj, the comparison between this and 
the spacelike method developed here follows, where we analyse their similarities and 
differences in applicability and efficiency. Finally, the two methods are combined in 
Section [5l to produce a more generally applicable method (as illustrated with the 
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recovery of the pertinent information about a "star" in AdS^) and we go on to look 
at extensions of the method to less symmetric cases in section [6l We conclude in 
Section [7] with a discussion and summary of the results. 



2 Background 

Recall the metric for AdS^ in coordinates {t,r,(p): 



ds^ = -f{r)de + ^^+r^d<l? (1) 



fir) = ^ + ^ (2) 

where R is the AdS radius. The existence of Killing vectors d/dt and d/d(j) leads to 
two conserved quantities (energy (E) and angular momentum (jjll)) and allows the 
geodesic equations to be written in the simple form: 

r^ + Veff = Q (3) 

where ' = ^ for some afRne parameter A, and Veff is an effective potential for the 
geodesies, defined by: 

K// = -(/W- + ^^-^) (4) 

where k = +1,-1,0 for spacelike, timelike and null geodesies respectively. Note 
that only null and spacelike geodesies can reach the boundary at r = cxd in finite 
coordinate time, and so these are the geodesies we work with when relating bulk 
physics to the boundary. The paths of a sample of null and spacelike geodesies 
through AdS-i are shown in figure [H where one observes that the null geodesies 
all terminate at the antipodal point on the boundarjU. This is in contrast to the 
spacelike geodesic endpoints, where there is a both an angular and temporal spread 
in their distribution, obtained by varying J and E (except in the S = case, which 
we focus on here, where the geodesies are all contained in a constant time slice). 
Consider a deformatiorH to the pure AdS spacetime by replacing ([2]) with: 

2 

fir) = l + ^-p{r) (5) 

where p(r) is an analytic function which is of comparable magnitude to at small r 
and tends to zero at large r. Now, in [lOj, the metric information was extracted by 



^Note that in .lOJ, the geodesic angular momentum was denoted L; here we use J to avoid 
confusion with Lt, which denotes the length of the system in the CFT (see section [3|. 

^This will not be the case in spacetimes which deviate from pure AdS-i, see figure [TT] in section 

m 

■^This is not the most general modification one could consider, however, in the more general case, 
one needs both null and spacelike probes to determine the metric, see section O 
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Fig. 1: A sample of geodesic paths in AdSs (with R — 1), all beginning at the same point 
on the boundary, with varying J and E. The null geodesies (left plot) all terminate at 
the same (antipodal) point, whereas this is not the case for spacelike geodesies (right 
plot). 
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Fig. 2: A static spacelike geodesic in AdS^ (left plot), with the regions A and B highlighted 
(right plot). 



using tiie endpoints of null geodesies and their relation to correlation functions in 
the field theory. Here we propose to use the endpoints of static spacelike geodesies in 
three dimensions, and the relation between their proper length and the entanglement 
entropy of a two dimensional CFT proposed in [5] to extract the bulk information. 



2.1 Entanglement entropy 

In Ryu and Takayanagi propose that the entanglement entropy Sa (in a CFT^+i) 
of subsystem A with (d — l)-dimensional boundary dA is given by the area law: 

Area of 

where ja is the static minimal surface whose boundary is given by dA, and gSJ+^^ 
is the Newton constant in (d + 2) dimensions. In the d = 1 case, 7a will be given by 
a geodesic line, and thus if we consider AdS^ with a (l+l)-dimensional CFT living 
on its boundary, and define two regions A and B on the boundary as in figure [21 Ryu 
and Takayanagi's proposal relates the proper length of the static spacelike geodesic 
shown to the entanglement entropy 5"^. Thus by considering a complete set of these 
geodesies, we can probe the entire spacetime from out near the boundary down to 
the centre at r = C^, as we discuss in the following section. 



^This assumes we are working in a non-singular spacetime; for the case where the central dis- 
turbance p{r) corresponds to that for a black hole, one can probe down to the horizon radius, rh, 
see section 
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3 Method for reconstructing /(r) 



Focussing on spacelike geodesies, and specifically those with zero energy (i.e. static), 
we have that: 



J' 



^ -/W(i-72 ) =0 (7) 



which can be combined with the angular momentum conservation equation J = r'^c, 
to give: 



This can then be re-cast as an integral equation along the geodesic path, where we 
note that the final angular separation will be a function of J only: 

Z<l>end fTmax 1 
d<P = 2 r^=d^ (9) 
r'start Jrmin f ^Jjjj^ ^ JL^ — \ 

where rmin is minimum radius obtained by the geodesic, and in the zero energy case 
is given simply by rmin = J ■ As the metric is divergent at the boundary r = oo, 
we introduce a cut-off T^ax and restrict ourselves to the region r <i r max 

n We also 

have that the proper length of the geodesic (also dependent only on J) is given by: 

C{J) = 2 dr (10) 



These two equations, ^ and (jlOp . will form the basis for our method of extract- 
ing the metric function /(r) at each r. 

Now, given that the spacetime in which we are working is asymptotically AdS, 
we can say that for r > r„ for some Tji which can be arbitrarily large (but still 
below the cut-off rmax), f{f) -|- 1 (with R set to one). Thus all static spacelike 
geodesies with angular momentum J > Jn = fn will remain sufficiently far from the 
central deformation p[r) such that they remain undisturbed by its effects, and in 
the limiting case J = r„ we can write: 

0n =2/ ,^^^ dr (11) 



vr / 2rl + {rl-l)r^ 



,1 n ' V / max i 

arctan = ) (12) 

2 \ 2r„Vr4„^ - (r2 - l)r^, 



max 



vr frl-l 



- - arctan ( — ) for Vmax > Tn (13) 



^This cut-ofF corresponds to the ratio between the UV cutoff (or equivalently the lattice spacing) 
in the CFT and the total length of the system: Vmax ~ Lr/a 
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where (j)n = 4>end — (f'start, and is the length of section B of the boundary in figure 
[2j Hence from the (p endpoints, which are specified by the our choice of region A in 
the CFT, we can determine r„ and we have that f{rn) = + 1. This wih be the 
starting point for an iterative method which will recover the metric from r„ down 
to zero (in the non-singular case). 

The naive way in which to now proceed is by taking a slightly smaller choice of 
minimum radius, r„_i < r„, and splitting up the relevant integrals in ^ and (|1U|) 
into two pieces, one from r„_i to r„ and one from r„ to Vmax- These integrals could 
then both be well approximated, the first by taking a series expansion about the 
minimum radius r„_i, and the second by approximating the spacetime as pure AdS, 
as in (jlip. We would thus end up with two simultaneous equations which could be 
solved to give r^-i and /(r^-i), and could then proceed in a similar fashion to 
obtain the the entire bulk metric, to an arbitrary level of accuracy determined by 
our choice of step size in r (which is determined by our choice of boundary region 
4'end — 4^ start)- However, it turns out there is a significant problem with this setup 
which prevents it being applied in practice. Specifically, the iterative process is 
unstable, with any errors in the estimates for r„_j and /(r„_j) leading to greater 
errors at the next step. This results in a rapid divergence of the estimate from 
the actual metric, and the iteration quickly breaks down. Whilst improving the 
approximations to the various terms in the integral can slightly improve matters, 
there is a way of avoiding this unstable setup (where we solve for the two unknowns 
simultaneously at each step) entirely, as we shall now demonstrate. 



3.1 Determining the angular momentum 

Consider the equations ^ and (fTOjl above; they both have very similar forms, and 
there is in fact a strikingly simple yet powerful relation between the two quantities, 
C and 4>- Taking the derivative of both with respect to J, the angular momentum, 
we have that: 



^=2 
dJ 



J 



- r2v7W(l-^) 



3/2 



dr 



and 



4 



drr. 



dJ 



(14) 



dJ 



J2\ 



3/2 



dr 



2 J 



4 



drr, 



dJ 



(15) 



Using the fact that J 



and noting that the divergent part of the integral 



cancels with the divergent second term in each equatiorO, we can see that the two 
equations are identical upto a factor of J, and we therefore have that: 



®It is straightforward to show this, and an equivalent calculation is given explicitly in the second 
appendix of [TUj . 
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dC d4> 
dJ^ dJ 



(16) 



which can be rewritten as 




mm 



(17) 



Thus we have the remarkable fact that the minimum radiuqil of the static space- 
Uke connecting any two points on the boundary is immediately calculable from the 
gradient of a plot of the proper length, C versus angular separation (j), see figure El 
This immediately provides us with one of the two unknowns we need at each step, 
and leaves us with only needing to calculate f{rmm)- This can be done iteratively, 
beginning at large r, by splitting up © (or ([TU]) ) and taking various approximations 
to each part of the integral, the full details of which are given in Appendix A. Unlike 
the original proposal for the method, this is very robust to any errors, and provides 
an efficient way of determining the bulk structure, as we see in the examples in the 
following section. 

The relation (I17p also allows us to more specifically determine the point at which 
the metric deviates from pure AdS; recall that on the first step of the iteration (with 
i = 0), we took the metric to be pure AdS, and after determining r„ using ([TT]) . set 
f{fn) = + 1) where we originally stated that r„ could be taken arbitrarily large. 
We can now explicitly check the radii at which the pure AdS assumption holds, as 
we can now determine the value of r^m corresponding to each (p separation of the 
endpoints, and hence plot 

ffi—i vs (f>n—i foi' each i. In pure AdS, we know that the 
relation is given analytically by r^m = cot ( ''^''"''"2'^°'°''* ), and at small enough angular 
separation, the two plots should coincide (this is also of course true on the plot of C 
vs (p, see figure [3|). This allows one to avoid beginning the iteration at an excessively 
large radius, which would reduce the efficiency of the extraction. 

We now address the issue of how confident one can be that the extracted solution 
matches the actual metric, before going on to consider some examples. 

3.2 Validating the extracted solution 

A natural question to ask at this point is on the uniqueness of the solution, i.e. 
is there more than one possible /(r) which gives the same boundary data for the 
geodesies? Then if there is a unique /(r), does this proposal for reconstructing the 
metric always find it, and not some alternative set of points (r„_j, /(r„_j)) which 
also solve equations (|50|) and (|56p without being the actual metric function? 

Considering the second question, it is quite simple to show that if the metric func- 
tion /(r) corresponding to the boundary data is unique, then the iterative method 

''Note that equation H17|) holds in any static, spherically symmetric spacetime; in those with 
less symmetry, such as angular variation of the metric as well as radial, one finds that the gradient 
^ gives the final angular momentum of the geodesic, but as this will not be conserved, it is not 
necessarily equal to rmin- 
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L[J] 




. 3: A plot of the proper length, £, vs the angular separation of the endpoints, 0, for static 
spacelike geodesies in an asymptotically AdS spacetime (red, lower curve), and in pure 
AdS (black, upper curve). The gradient, dC/dcj) at each point provides the angular 
momentum, J, for the corresponding geodesic. When the angular separation is small, 
the geodesies remain far from the centre, away from the deformation, and hence both 
curves coincide. 
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must recover it (up to a level of accuracy determined by the number of steps). We 
will show that if this is not the case, then either the metric function was not unique, 
contradicting our assumption, or the estimate does not in fact correspond properly 
to the boundary data. 

Take the extracted points {vn-i, f{rn-i)) for i = 0, . . . , n, and use them to con- 
struct an interpolation function, which is then our estimate for the metric function. 
We can then use this estimate to compute the proper length and angular separa- 
tion of all spacelike geodesies passing through the spacetime. If the generated data 
matches with the original data from the field theory, we have successfully produced 
an estimate for an actual bulk metric, and by our assumption of uniqueness, this 
function must be f{r). 

If the generated data fails to match correctly to that from the field theory, we 
can deduce that we haven't in fact produced an estimate for /(r), but instead that 
our (r„_j, /(r„_j)) are simply a set of points which solve the equations (f50]) and 
(j56l) . In this case, the iterative step size used to produce the estimate was too large, 
and the extraction procedure should be repeated with a smaller step size. Once 
the new estimate has been produced, the above test can again be applied; this can 
continue until an actual estimate of /(r) is recovered. 

Finally, one should note that at an infinitesimally small step size, one will use the 
completqj set of spacelike geodesies to probe the spacetime, generating a continuous 
estimate for /(r) from r„ down to zero. As such the data generated from our 
estimate must correspond to that from the field theory, as it was all used in its 
production. Thus, by uniqueness, the estimate must correspond to f{r). 

A basic argument for the uniqueness of the bulk metric corresponding to the 
field theory data (in our case, the proper length of the static spacelike geodesies as a 
function of the angular separation of the endpoints) follows from a comparison of the 
local degrees of freedom on each side, by noting that this data and the geometry of 
the constant time slice we wish to recreate contain the same amount of information, 
as /(r) is a function of the radial coordinate only. When coupled with the knowledge 
asymptotic behaviour of the spacetime (that it approaches pure AdS at large r), we 
have the boundary conditions needed to ensure that the metric function is unique. 
In less symmetric cases one has more freedom in the metric, but correspondingly 
one also has more information with which to determine this, see section [6] for further 
comments on these scenarios. 

Finally, one should note that this is simple argument does not constitute in 
any way a proof of the existence or uniqueness of the solution, as here the focus is 
on demonstrating how an intriguing link between field theory and the bulk leads 
to a remarkably simple process for calculating numerically the corresponding bulk 
metric. With this in mind, having argued that with suitable checks the extracted 
solution should be an estimate for /(r), we now proceed to some examples where 
we examine the accuracy of such estimates. 

*By complete, we mean all geodesies which have minimum radius r,„i,i < r,i, where r,i can be 
taken arbitrarily large 
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step size a (4) (5 (1) 7 (8) X (3) i] (2) A (1) 



0.1 
0.05 
0.01 
0.005 



3.75 
3.81 
3.94 
3.95 



0.70 
0.79 
0.85 
0.93 



7.99 
7.95 
8.19 
8.01 



3.03 
3.02 
3.01 
3.01 



1.99 
1.99 
2.00 
2.00 



1.00 
1.00 
1.00 
1.00 



Table 1: Best fit values (to 2 d.p.) for the /fiti(r) parameters a, (3, 7, rj and A, with the 
actual values indicated in brackets. 



3.3 Examples 

To illustrate the procedure for metric extraction, we begin by considering some 
examples of deformations of the pure AdS metric. In the cases considered we have 
taken the proper length and angular separation of the endpoints to be known from 
the relevant field theory, and taken a linear step size in J (and. hence T^nin 

). The 

method of Appendix A is then applied for a variety of step sizes, and the resulting 
estimates for /(r) are plotted alongside the actual curve. The three deviations from 
pure AdS we consider are the following: 

„ , , 9 4r^ 3rsin(2r) , , 



10sin2(3r) 
r3 + 1 



/2(r) = l+r^ + (19) 



. / N 9 10sin^(10r, , , 

/gr =l + r' + 3XT^ 20 

r'^ + 1 



where each gives a non-singular, asymptotically AdS spacetime. These functions 
were chosen as tests of the extraction method because they provide clearly visible 
deviation from the pure AdS metric of /(r) = + 1. The first example also 
corresponds to one used in [10] in an alternative method for extracting the bulk 
information (see section H]), and despite the similarities between /2(?') and fsir), we 
shall see a noticeable difference in the accuracy of their extraction for larger step 
sizes. 

For the first example we use four choices of step size in r, namely Ar ~ 
0.1,0.05,0.01 and 0.005, and compare the accuracy of the generated curves to the 
actual function; this is done by considering best fits to the numerical estimates, 
obtained by using a non-linear fit to the following function: 

r r ^ 1,2 , Xr sin(r?r) 

to give values for the various parameters. The results are shown in Table [U with 
the corresponding data points plotted in figures H] and O 
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Fig. 4: The data points for the largest two step size estimates for /i(r), compared with the 
actual curve (in blue). Whilst both give good estimates to the curve, the step size of 
0.1 (left) deviates at a higher r than when using a step size of 0.05 (right). 




Fig. 5: The data points for the next-to-smallest step size estimate for /i(r), compared with the 
actual curve (in blue). The fit here appears very good even close to r = 0, however, 
Table [1] shows that we still need to go to a smaller step size in order to accurately 
extract values for a, (3 and 7. 
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^ . . . . , . . . . , . . . . : . . . . ^ 

1 2 3 4 5 



^ — . — . — . — . — , — . — . — . — . — , — . — . — . — . — , — . — . — . — . — ^ j- 
1 2 3 4 5 



Fig. 6: The data points for the largest two step size estimates for f2{r), compared with the 
actual curve (in blue). Despite the larger deviation from pure AdS than in example 1, 
both the estimates here provide good fits to the curve. 



Step size 


X (10) 


r?(3) 


A(l) 


0.1 


10.32 


2.99 


1.06 


0.05 


10.08 


3.00 


1.01 


0.01 


10.05 


3.00 


1.01 



Table 2: Best fit values (to 2 d.p.) for the /fit2(7') parameters x, rj and A, with the actual 
values indicated in brackets. 

Prom Table dl which contains the data for the estimates of /i(r) we see that 
there is a very good fit to the actual values of the six parameters from our non- 
linear fit (|2ip . even at the largest step size we consider. Indeed, by eye it is hard to 
tell any difference between the accuracy of the estimates except at very small radii. 
This is mainly due to the relatively slow variation of /i(r) with r, which ensures the 
various approximations we take in order to produce the estimates remain good even 
for the larger step sizes. Whilst it appears that taking a smaller step size is rather 
superfluous, it should be noted that the finer structure parameters (namely a, j3 and 
7) would need the smaller step size data in order to be determined to a high level 
of confidence. Our choice of non-linear fit function is also rather specifically chosen 
to match the example; if one did not know beforehand the form of fi{r) one would 
want to take smaller step size estimates in order to obtain data down as close to 
r = as possible (as is discussed at the end of the section), to ensure that any finer 
structure was not being masked, and also as a check on the validity of the previous 
estimate. 

We see similar behaviour in the second example, where we have chosen a slightly 
more fluctuating function to attempt to recover. Here we use the three largest 
choices of step size in r, and the data generated in each estimate is shown in figures 
[6] and [71 where we also include a plot of the actual function f2{T) as comparison. 

We can again use a non-linear fit to evaluate the estimate; in this case we use a 
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1 2 3 4 5 



Fig. 7: At a step size of 0.01, the estimate data for /2(f) matches the actual curve (in blue) 
almost exactly, even close to r = 0. 



Step size 


X (10) 


V (10) 


A(l) 


0.1 


7.49 


8.03 


0.29 


0.05 


11.60 


10.00 


1.25 


0.01 


9.96 


9.99 


0.99 



Table 3: Best fit values (to 2 d.p.) for the ffit2{r) parameters x, V and A, with the actual 
values indicated in brackets. 

function of the form: 

/,,M = l + r= + ^^=5!^ (22) 

and the results are shown in Table [5J 

Thus far everything is progressing as expected: the smaller step sizes are pro- 
ducing closer fits to the curve, and better estimates for the values of the various 
parameters. In these first two examples, we even have that the largest step sizes 
produce good fits to the curves; do we ever see a large increase in accuracy over 
our choice of step size? If we consider the third example (which was obtained by 
increasing the value of 77 from the second example), where the function oscillates 
more wildly at low r, we do see a significant improvement in the estimates as the 
step size decreases. Proceeding as before, we see that for the largest step size of 0.1, 
the method has difficultly in following the rapid oscillations at low r; this is then 
significantly improved upon in the subsequent estimates, as shown in figures [8] and 
El and in the non-linear fit data given in Table O 

As expected, the smaller step size again produces a closer fit to the actual curve, 
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Fig. 8: The data points for the largest two step size estimates for fsir), compared with the 
actual curve (in blue). The reduction in step size from 0.1 (left) to 0.05 (right) gives 
a marked improvement in the fit of the points to the curve at low r. 




Fig. 9: The data points for the smallest step size estimate for f^{r), compared with the actual 
curve (in blue). This level of precision gives a very good fit to the curve, and this is 
mirrored in the highly accurate estimates for the function parameters, given in Table[3] 
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however, in this third example, the largest step size fail to give accurate data for 
the unknowns x, rj and A, although it does make a reasonably close fit to the curve 
until the iterative process breaks down. 

Finally, one should comment on the fact that the deviation of the estimate from 
the actual curve does not apparently prevent the iteration from continuing to give 
sensible looking (although erroneous) values in subsequent steps. Whilst appearing 
to allow for an incorrect determination of the metric, applying the checks described 
in section 13.21 (reconstructing the field theory data using the metric estimate) will 
quickly highlight any areas in which the estimate for /(r) has deviated from the 
correct function. As stated before, this merely indicates that the step size in r 
was too great for the iterative method to properly be effective in extracting the 
information using the approximations chosen in Appendix A. Aside from simply 
reducing the step size, or using better approximations (such as at each step creating 
an interpolating function estimate for /(r) using the already determined data), there 
are other possible resolutions of this problem to further optimise the extraction. One 
could take either a non- linear step size in r to include more terms near r = 0, or 
simply take appropriately varying step sizes depending on the fluctuations of the 
metric; where the metric is varying rapidly with r the step size could be reduced. 
Thus by making several passes, reducing the step sizes at appropriate r each time, 
the estimate of /(r) can be significantly improved without considerably increasing 
the computation time. 

We now conclude the examples section by briefly investigating how the method 
is affected in spacetimes with a wildly fluctuating interior, and how one can apply 
the above to maintain a high degree of accuracy. 

3.4 Maintaining accuracy in wildly fluctuating spacetimes 

The third example of the previous section has shown that in wildly fluctuating 
spacetimes one needs smaller step sizes in order to guarantee accuracy of the estimate 
for /(r) down to small r. Here we provide two further examples to show how the 
method breaks down if the frequency of the fluctuations is sufficiently increased, 
and how one can adjust the step size to compensate. 

Firstly, one observes that it is not simply the frequency of the oscillation which 
causes the extraction to break down, but also the amplitude; this can be seen in 
figure [TOl where the estimate continues to follows the actual curve closely whilst the 
amplitude of the oscillations is small. The two examples shown in the figure come 
from considering modifications to example 3 where the sin^(lOr) term is replaced 
by first sin^(20r) and then sin^(30r); as stated, one still obtains a relatively good 
fit to the curve using the smallest step size, although in the more rapidly oscillating 
case the fit does deviate slightly more from the correct curve, especially near the 
peaks at low r. 

This behaviour is important, as it means that even in metrics with a large and 
rapidly varying interior, one can use a reasonable step size to extract the metric with 
confidence down to a fairly close distance to the centre. After checking the estimate 
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Fig. 10: Plots of /3(r) with the sin^(lOr) term replaced by sin^(20r) (left) and sin^(30r) 
(right), along with estimates generated with a step size of 0.01. Interestingly, whilst 
sufficiently increasing the frequency of the metric oscillations does reduce the depth 
to which the metric is accurately extracted, it does not adversely affect the accuracy 
of the fit to that point. 



by recreating the field theory data, one can then continue the extraction from that 
point with better approximations, and a smaller step size (beginning slightly further 
out than the final terms so as to give some overlap with the initial estimate and check 
the consistency of the estimates) in order to fully reconstruct the metric function. 

In any case, the more exotic spacetimes one might wish to consider may not have 
only one independent metric function /(r) to extract, and in order to fully determine 
the metric in these more general cases, one may also need to consider the use of null 
geodesic probes. Thus having now established the principles of the method, and 
seen some examples, we go on to look at comparisons with an alternative method 
of extracting the bulk metric proposed in previous work. 

4 Comparison with an alternative approach to metric 
extraction 

After seeing in the previous section examples of how the extraction works in practice, 
we now consider how this method (S) based on spacelike geodesies compares to an 
alternative method involving null geodesies (N). Before we do so, however, we firstly 
provide a short review of this different approach to probing the bulk, which was 
originally presented in |10j . 

4.1 Review of the null geodesic extraction method 

For a spacetime of the form of ([T]) with metric function /(r) as in (|18|) say, we 
can consider the full set of null geodesic paths through the bulk, which is obtained 
by choosing some arbitrary starting point on the boundary and varying the ratio, 
y = J/E from zero to one, see figure [TT] 

From this plot of the endpoints, if one takes the gradient dt/dcj) at any point, one 
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Fig. 11: Null geodesic paths passing through a modified AdS spacetime, all starting from the 
arbitrary point t — 0, (p = on the boundary and with y > 0. The corresponding 
full spectrum of null geodesic endpoints for this spacetime is shown on the right. 
(Redrawn figure from [10| ) 

obtains the value of y for the corresponding geodesic. This is in a sense the "extra" 
piece of information (analogous to considering dC/d(p, see section \3A\i determined 
from the CFT which allows the geodesic probes to extract the bulk metric; here it 
is the ratio y of angular momentum to energy which is obtained, in the spacelike 
geodesic method it was simply J. After determining the first term of the iteration 
by taking the spacetime to be pure AdS far away from the centre, one can then 
take similar approximations to those given in Appendix A to split up the relevant 
geodesic equation: 




and combine with the relation y = rmin/ \/ f{fmin) to iteratively extract the met- 
ric. At this point it is worth making a computational observation about the two 
approaches; both involve almost identical procedures for iteratively extracting the 
metric, and as such are of comparable efficiency. There are, however, a number of 
fundamental differences between them, as we shall now discuss. 

4.2 Dimensional applicability 

Whilst in method N was applied to the specific case of AdS^, it is equally 
applicable in an arbitrary dimensional spacetime, AdSn+i (for n > 2), assuming 
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one could obtain the endpoint information from the appropriate field theory on 
the boundary. Whilst the principles of method S can also be applied in arbitrary 
dimensions, it is no longer clear as to whether the proper length of the spacelike 
geodesic is so readily extractable from the CFT in anything other than the re = 2 
case. In higher dimensions, the area of the minimal surface which corresponds to 
the entanglement entropy is no longer the length of a spacelike geodesic, and the 
method would need to be modified to take this into account. This could be achieved 
either by using some expression for the minimal surface instead of the proper length 
equation ([TU]) . or by demonstrating an alternative route to determining the proper 
length. 

4.3 Singular spacetimes and those with significant deviation from 
pure AdS 

One of the main limitations of method N is that it cannot probe past a local maxi- 
mum in the effective potential for the null geodesies (see figure [12]); it cannot there- 
fore probe close to the horizon of a black hole for instance. The method presented 
here would have no such problem, as the spacelike geodesies can reach arbitrarily 
close to the horizon while still being able to return to the boundary. For example, 
consider a five dimensional Schwarzschild-AdS spacetime with metric function /(r) 
given by: 



where we have set = R = 1. As was shown in [TO], using method N one is only 
able to probe down to a radius of r = 2, as at this point the effective potential 
for the null geodesies becomes a local maximum. Method S, however, allows the 
bulk information to be fully determined all the way to the horizon radius, = 1. 
Similarly, for those non-singular spacetimes with large enough deviation from pure 
AdS so as to allow for null geodesic orbits (the signature of a non-monotonic effective 
potential), one has no problem extracting the full metric using method S, as in the 
second and third examples of section 13.31 

4.4 The overall conformal factor 

Finally, one should point out that the method presented here is also sensitive to the 
overall conformal factor of the metric, whereas method N is not. This simply stems 
from the fact that for null geodesies, ds"^ is zero, and hence for any metric: 



the conformal factor immediately drops out. For spacelike geodesies however, ds"^ = 
1, and thus the Q{r) term remains. Whilst this conformal factor ^}{r) presents us 
with another unknown, we shall see in the following section how it can be determined 
by combining the two methods (N and S) together. 



/(r) = l+r2-4 



(24) 



r 




(25) 
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Fig. 12: Plot of the effective potential for two null geodesies with similar y, in some arbitrary 
spacetime. The null probe which follows the solid effective potential will go into 
circular orbit due to the local maximum; the geodesic with slightly lower y (dashed 
red line) then has significantly lower rmin, and this finite jump in the minimum radius 
causes the iterative extraction method to break down. 

5 Applying the two methods together 

Having compared the relative merits of the two methods, we now investigate how 
it is possible to use them in conjunction with one another to determine the metric 
in more general cases. Thus far we have restricted ourselves to considering metrics 
of the form of ([I]), however, we can look to extend this further by considering the 
most general static, spherically symmetric spacetimes, given by metrics of the form: 



-fir)df + h{r)dr^ + 



(26) 



where we have incorporated the conformal factor ri(r) of ([25]) into two new functions 
/(r) and h{r) (and rescaled the radial coordinate accordingly). Using either method 
independently to recover the metric fails because of the presence of three unknowns: 
r, f{r), and h{r) with only two independent equations with which to determine 
them. We can, however, use both methods in conjunction, as outlined below, where 
we restrict ourselves to the (2+l)-dimensional case in accordance with section ^^21 
For a spacetime of the form of (j26|) , we have the two constraints on the energy 
and angular momentum from before: 



E = f{r)i 



(27) 
(28) 
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along with the modified expression involving the effective potential: 



.2 / E'^ J 



\h{r) ^ f{r)h{r) h(r)r'^ ) ° ^^^^ 

We immediately see that for the zero energy spacelike geodesic paths we do not 
obtain any information about the function /(r) (as we would expect, as /(r) affects 
the time coordinate, which is kept constant in the E = case), and our integrals 
for the separation of the endpoints and proper length are given by: 



r 



"Jend - <Pstart = 2 / , dr (30) 

" rJC-1 



and 



£ = 2/ dr (31) 



rmn i / 1 — ^ 



We can thus use the static spacelike geodesies to determine h{r), from r = 
to an arbitrarily large r„, by applying the extraction method proposed in section [3] 
and Appendix A. Specifically, for each we have the corresponding h{ri), and from 
this one can generate a best fit curve, hf[t{r). One then is left with extracting the 
/(r) information from the null geodesic endpoints: for a null geodesic in a bulk with 
metric (pUj) . we have that 



'■'\t = 2r —^^=,r (32) 



with the minimum radius given hy y = rmin/ \J f{fmin)- If we now replace the 
function h{r) with our estimate /ifit(r), this becomes 



-y^) ^..^ (33) 



tstart Jrrmn f (r) . 



y 

f{r) 

which contains only two unknowns, as the parameter y is given by the gradient of 
the endpoints (see section HJ]) . We can then use the iterative method of [lOj (the 
relevant equations are given in Appendix C) to recover the second metric function, 
/(r), and the bulk information has been extracted, as we see for the two examples 
which follow. The main area of concern would be whether significant errors in 
recovering /(r) appear unless the estimate function for h{r) is highly accurate; one 
can investigate whether this is so by running the extraction of /(r) several times 
using a different estimate for h{r) each time. We see how this affects the results 
in the first example below. Finally, one should note that the depth to which the 
metric can be recovered is subject to the same restrictions as before (see section 
T3l): for example in singular spacetimes, whilst the spacelike geodesies can probe 
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down to the horizon radius, (and we thus obtain h{r) down to h{rh)), the null 
geodesies can only probe as far as the first local maximum in the effective potential 
(figure [T2]) ■ at some > rh, leaving /(r) undetermined for r < r/j2. Nevertheless, 
by combining the two different approaches to probing the bulk, we have obtained 
more information than is possible using either individually. 

5.1 Example 1: Testing the combined extraction procedure 

Consider a spacetime where the two metric functions /(r) and h{r) are given by the 
following: 



„, , , 9 7r 2rsin(5r) , ^, 

fir) = 1 + r2 - + \ J 34 



, , , / o 4r2 3rsin(2r)\ , , 

h(r) = 1 + - — — + — — L (35) 

^ ^ V (r2 + l)(r2+8) r^ + l y ^ ^ 

Whilst this is in no way meant to be a representation of any physical deformation 
of the bulk, it is a good test of the combined extraction method, as it provides a 
monotonic effective potential for the null geodesies, and so allows us to probe down 
to r = 0. One can also use the similarity between this spacetime and that described 
in the first example of section [331 namely that we have h(r) = fi{r)~^. This was 
deliberately chosen so the part of the metric probed by the spacelike geodesies is 
exactly as it was in the case of example 1; the change in /(r) has no effect on the 
results, and thus the best fit estimates for h{r) are exactly those specified by the 
values of the parameters in Table [TJ We therefore have four different estimates for 
h{r) (one for each of the four choices of step size used), and we label them hQ i{r) 
through to hQQQ5{r), where the subscript refers to the step size. All that is left to 
do is to attempt to recover /(r) via the null geodesic dat4j for each fit to h{r), and 
compare it firstly to the actual values, and also to those obtained using the exact 
function h{r) rather than an estimate. The results are analyzed using a best fit of 
the form of ()21|) and are presented in Table HI 

We see quite clearly from the table of results that even using our roughest esti- 
mate for h{r), namely hQ i{r), we obtain a highly accurate estimate for f{r). Indeed, 
the limiting factor is not the accuracy of the estimate for h{r), rather it is the choice 
of step size and starting y in the null geodesic part of the extraction (see footnote 
i). 

5.2 Example 2: Radiation in AdS^, a toy model 

As the two extraction methods give such good fits when applied sequentially, we 
now turn our attention to a less trivial example, where we consider a gas of radi- 



®As we saw in [TD], one can use a range of different step sizes in y to obtain varied levels of 
accuracy in the metric extraction; as we are not intending to specifically analyze the null geodesic 
method here, we simply choose a starting value of y = 0.9985, and a step size of Ay — 0.0005, as 
these are sensible values for the example given. 
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a (7) 


P (1) 


7 (13) 


X (2) 


T] (5) 


A (15) 


ho.iir) 


6.81 


1.03 


12.49 


2.00 


4.99 


14.92 


ho.obir) 


6.81 


1.03 


12.48 


2.00 


4.99 


14.92 


ho.oi{r) 


6.80 


1.03 


12.48 


2.00 


4.99 


14.92 




6.80 


1.03 


12.48 


2.00 


4.99 


14.92 


h{r) 


6.80 


1.03 


12.48 


2.00 


4.99 


14.92 



Table 4: Best fit values (to 2 d.p.) for the ffit{r) parameters a, (3, 7, x. V and A, with the 
actual values indicated in brackets. We see that even our roughest estimate for h{r) 
is close enough for the extraction of /(r) to be highly accurate. 



ation in AdS^. There have been numerous papers exploring this and other closely 
related geometries in various dimensions, such as [HI HSJ [TSj Ej , and we focus here 
purely on our ability to recover the metric information via our numerical extraction 
methods. Firstly, we note that whilst restricting ourselves to three bulk dimensions 
does make our spacelike geodesic method fully applicable (see section B?^ . it also re- 
stricts the physical realism of the model due to the non-dynamical nature of gravity. 
Nevertheless, it provides a good toy model for radiating "stars" in AdS spacetimes, 
and allows us to demonstrate how well the pertinent information (e.g. the "star's" 
mass and density profiles) about the bulk can be recovered. We consider a perfect 
fluid solution to Einstein's equations, with the pressure P{r) set equal to half the 
density, p{r)/2, as for radiating matter the stress-energy tensor is traceless. For a 
metric of the form of (1261). we find thatrj: 



h{r) = (1 -m(r)) ^ (36) 

and 

/«=(fy) (37) 

where the mass function is defined by: 

m{r) = 2 / p{r)rdr (38) 







and Poo is the leading coefficient of p{r) at large r, and is given by poo ~ p{r)r^ as 
r CO. We obtain from the field equations a pair of coupled ODEs for m[r) and 
p{r): 

m'{r) = 2p{r)r (39) 
6 + 3,(r) _^VM^^ 



1 + — m(r) p[r)r 



^"We set R = 1 and SttGs = 1 for convenience. 
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Fig. 13: The density and mass profiles (top plots) for a "star" with central density po = 0.8, 
along with plots of the corresponding metric functions /(r) and h{r) (bottom). 

which when combined with the relevant boundary conditions m(0) = and p{0) = po 
can be numerically solved to allow us to generate the geometry of the spacetime (see 
figure [T3|) . The condition p(0) = po specifies the internal density of the gas, and po 
is the single free parameter of the system: pure AdS is recovered when pQ = 0. 

Before we begin with the metric extraction, we should make a comment about 
the features of such spacetime at large radius, as there are significant differences in 
the asymptotic behaviour of the metric depending on the choice of pQ. For pQ ^ 0, 
we have that the asymptotic behaviour of the metric functions is given by 

h{r) ^ {l + r"^ - M)~^ and /(r) ^ 1 + - M as r ^ oo (41) 

where M > is a constant. If M > 1 we have that the metric becomes the BTZ 
black hole solution at large r (see [15 1 116 1 ^oi example] for more details); this poses 
a problem for the method involving null geodesies, as we can no longer use them to 
probe the full range of r. Whilst this is due to the form of the effective potential 
(see figure [T^ , it is not due to the local maximum problem we saw in section 14.31 
Rather here we no longer have geodesies which can usefully probe the spacetime 
away from the centre: for the full set of null geodesies (obtained by varying y for 
zero to one), the minimum radius reached by the geodesies is bounded from above. 
We thus cannot take r^m to be arbitrarily large on the first step of our iteration, 
which was necessary for us to begin extracting the metric (although we should note 
that we could still apply the spacelike method to extract h{r) in this scenario). 
Instead however, we will consider the region < M < 1, corresponding to conical 



24 




Fig. 14: Effective potentials for null geodesies in a spacetime with M = 8. The upper (red) 
potential is for y = J/E = 0.9999; no matter how close to one the ratio J/E 
becomes, the minimum radius (defined by Ve/f = 0) remains small. 



defects, in which both methods are apphcable and is obtained by taking po to be 
small O 

Let us then proceed with recovering the metric in the specific example shown in 
figure [T3l where we have set po = 0.8. Bearing in mind that our goal is to firstly 
reconstruct the functions /(r) and h{r), and then use these to determine the mass 
and density profiles (m(r) and p{r) respectively) of the star, we begin by applying 
the spacelike geodesic method (with step sizes of 0.1, 0.05 and 0.01) to produce three 
estimates for h{r), the most accurate of which, namely ho,oi{r), is shown in figure 
[TE\ Whilst in the previous example we defined /i(r) explicitly by hand, and so knew 
the form of the function with which to apply the non-linear fit to generate the best 
fit curve /ifit(r), here we do not have such a starting point. Instead, we use the data 
points (r„_j, /i(r„_j)) to generate an interpolating function which will serve as our 
h{it{r). Thus although we cannot write down an explicit form for hfii{r), we can use 
the interpolating function to then carry out the next part of the extraction process, 
namely using the null geodesic probes to recover f{r). 

Using the third (and most accurate) estimate for h{r) in the modified null 
geodesic method of section [5] and Appendix C, we produce the estimate for /(r), 

^^One should also note from (|41|l that our iterative equations for recovering the metric need to be 
modified to take into account the new asymptotic behaviour, as we no longer have that the metric 
is given by /(r) ~ + 1 at large r. Thus we say that for r > r„ we have that /(r) and h{r) are 
given by (I41|l . and modify the approximations to the integrals for (pn-i and £„-i accordingly. 
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Fig. 15: The third (and most accurate) estimate for h{r), where the fit is good down to 
r ~ 0.1 (left plot). The estimate for /(r) generating using this approximation to 
h{r) is given in the right plot, and we see that it too appears accurate down to very 
low r. 



f{it{i^), also shown in figure [TSi we have now reconstructed the star metric. Although 
if we so wished we could have taken smaller step sizes to improve both the estimate 
of h{r) and that of /(r), we now continue with the ones we have. 

How do we use the metric functions to determine the mass and density infor- 
mation for the star? From (|36p it is immediately obvious: we can rearrange the 
equation to solve for m(r), and substitute in our interpolating function h{it{r) to 
give an estimate for the mass profile: 

mfit(r) = l + r2--i— (42) 
hfitir) 

and we obtain a fit for the density profile in similar fashion, by using the above 
estimate for m(r) in ([5^ . to give: 

, , mL(r) , , 

m = 43 

These two fits are plotted against the actual functions m(r) and /9(r) in figure [T6l and 
we see that by using the metric function data /ifit(?') we have obtained reasonably 
good estimates of the mass and density profiles of the star, aside from at very small 
r, where the errors from the estimate of /i(r) become noticable. What is noticeable 
is that the estimate for p{r^ fails at higher r than any of the others; this is due to 
the use of the derivative of the interpolating function mfit(r) in its construction, and 
is dealt with later (see below). 

One now asks the obvious question of why it was necessary to extract the function 
/(r) at all, seeing as we have apparently just reconstructed the information about 
the star simply by using /ifit(?")- This is where we recall that we should be assuming 
that a 'priori we knew nothing about the origin of the metric's deviation from pure 
AdS. In fact, this has not been the case. Whilst our expressions for /i(r) in terms 
of m(r) and the mass m(r) in terms of the density /o(r), ()36p and (I39p . stem from 
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Fig. 16: Estimates for the mass and density profiles for our "star" . As with /lo.oi(^) and /fit(r), 
these match the actual curves closely until low r, although the density estimate Pfit(r) 
(dashed) fails at noticeably higher r than the others. Included in the lower plot are 
alternative estimates for for the density profile, obtained from (P7|) (closest fit) and 
61 (solid) 
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the dimensionality of the bulk (e.g. in higher dimensions one would have the m(r) 
term multiplied by some negative power of r), in defining /(r) by (j37p we have 
already taken the matter content to be a gas of radiation, which sets P{r) = p{r)/2 
and removes the pressure profile as an unknown. Given this knowledge, one could 
indeed have simply used the spacelike geodesic method say to extract the information 
about the star, as h{r) gives m{r), and m(r) gives p{r). Extracting /(r) becomes 
a necessity, however, if one removes the assumption about the matter content; then 
one also has to compute the pressure profile. It is most easily determined (once we 
have our fits for /(r) and h{r)) from the Grr component of Einstein's equations, 
and we have that: 

Pm{r) = ^ Jf^"^ , , -1 (44) 

which in our example corresponds to pfit{r)/2. Therefore by also plotting 2Pfit(r) 
in figure [T6l we can see how close the fits generated by the two different expressions 
(|43p and (j44p match, and this provides a check that the matter content is indeed 
that of a gas of radiation and confirms that our expression, ([37]), for /(r) is correct. 
Interestingly, we see that this expression provides a slightly better fit to p{r) at small 
r than that from ()43p . This is simply because ()44p includes ffnif) terms, and the 
non-linear step size in r in the null extraction method generates a greater amount 
of data points at low r for the estimate for f{r), thus allowing the derivative of the 
interpolation function to be more accurately determined. We can obtain the best fit 
at low r by using /fit(r) in ([57|) and solving for p{r) (see figure where we have 
avoided using derivatives 

Finally, we can use the estimates to give a numerical value for our free parameter 
Pq. Taking poo as having been calculated from the asymptotic fall off, and approx- 
imating the value of /(O) as 0.525, we obtain a value of 0.76, compared with the 
actual value of po = 0.8. Whilst the match is fairly good, this is where the accuracy 
of the estimates for /(r) and h(r) become very important; in taking /(O) = 0.525 
we have discarded the final few iterations of f{it{r) at small r, which lead to a kink 
in the curve, as being erroneous and due to an incomplete recovery of h{r). This is 
a reasonable assumption to make, as in our previous examples we saw that for too 
large a step size the method of generating hfit{r) fails to reach down to r = 0. We 
also have the data from the higher step size fits (/io.i('^) and /io.05(^)) with which to 
analyse the accuracy of our estimates for h{r) at low r. However, as it is the small r 
region from which the numerical value of po is calculated, in order for it to be con- 
fidently extracted one must ensure the estimates hf^tir) and ffit{r) are thoroughly 
checked for r close to zero. 



^^One should note this does firstly require the value of poo to be determined from the fall off of 
p(r) at large r; this is however available from our earlier fit to p given in (|43p . 
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6 Extensions to less symmetric cases 



In all of the above we have taken the spacetime metric to be both static and spheri- 
cally symmetric, however, we now consider how the methods for extracting the bulk 
presented here could be extended to include more general cases. 

Reducing the amount of symmetry removes conserved quantities from the geodesies; 
spherical symmetry gives us conservation of angular momentum, time translational 
symmetry gives us energy conservation. Consequently, there will be additional un- 
knowns introduced in our analysis of the geodesic path, as we will need to know more 
details about its route through the bulk; this should not prove a problem, however, 
as there will also be further information available from the geodesic equations. 

Consider for example the non-spherically symmetric (but still static) case. Be- 
fore, when there was no angular dependence in the metric, we considered a series of 
geodesies which probed deeper and deeper into the bulk - in other words, we had 
one which probed down to each r„_j. These were specified by the angular separation 
of the endpoints on the boundary, and the actual values of the (pstart and (pend were 
unimportant. This allowed us to reconstruct the bulk step by step, one value of 
f[rn~i) at a time. 

Now, what is the analogous method in the non-spherically symmetric case? At 
each step of the iteration we can still consider some fixed angular separation of 
the endpoints, however, we must also vary (j)start from to 27r (with some choice 
of slicing sufficient to give an accurate estimate), such that for each iterative step 
we recover a "ring" of information about the metric. The subsequent steps then 
recover smaller and smaller rings, extracting the metric function down to the centre 
of the spacetime. This is the basic idea of the extraction method; finalising a more 
detailed procedure which gives high accuracy within a reasonable computational 
time is subject of current research, one now has two step sizes to consider: the 
slicing in (f) and the radial steps in r. 

Finally, we should recall that higher dimensional cases offer further complica- 
tions, as mentioned in section WTH as although the null geodesic method is already 
applicable in such cases, the spacelike method is not. In principle though, the ideas 
still hold; one would need an expression for the correct minimal surface correspond- 
ing to the entanglement entropy (see [9] for more details on this point) which could 
then be treated in much the same way as the geodesic proper length, as they will 
each probe to a certain depth in the bulk, and those remaining at large r will be- 
have as in pure AdS. Completing the analysis for these cases is again the subject of 
further research. 

7 Discussion 

In this paper we have seen how the bulk geometry can be extracted (in asymptot- 
ically Anti-de Sitter spacetimes) using the entanglement entropy information ob- 
tained from the corresponding boundary field theory. In the three dimensional case. 
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the entanglement entropy of a subsection A of the 1 + 1 dimensional boundary is 
proportional to the proper length of the static spacelike geodesic connecting the 
endpoints of A (see figure [2]). Using this relation, Q, together with the endpoint 
data allows both the minimum radius, rmin, of the spacelike geodesic and the value 
of the metric function f{rmin) at this point to be determined, provided sufficient 
information about the bulk is known for r > r^ain- Thus by starting from large r, 
where the metric is approximately pure AdS, one can probe further and further into 
the bulk using geodesies connecting progressively smaller regions on the boundary. 

This extraction of the metric is made significantly more straightforward by an 
interesting relationship between the proper length of the geodesic and the angu- 
lar separation of its endpoints. Specifically, the gradient, dC/dcj), gives the angular 
momentum, J, of the corresponding geodesic, which in the static, spherically sym- 
metric cases considered here, is equal to the minimum radius the geodesic probes 
down to in the bulk. 

An iterative method for recovering the metric information in practice by approx- 
imating the relevant geodesic equation was thus then proposed, and a number of 
examples were given. The iterative method was developed in analogous way to the 
method presented in [TU], which used the endpoint data of null geodesies to similarly 
extract the bulk information, and was reviewed in section 14.11 

A comparison of the two methods was then made, investigating their relative 
advantages and disadvantages; this highlighted a number of differences in their rel- 
ative applicabilities. Whilst the method of |10], which uses null geodesies as probes, 
runs into problems when encountering geometries with significant deviation from 
pure AdS (as this leads to a non-monotonic effective potential for the geodesies 
which limits the depths to which the metric information can be recovered), no such 
limitations arise for the method given here involving spacelike geodesies, unless the 
metric is also singular. On the other hand, the null geodesic method is completely 
applicable in any number of dimensions, whereas although the principle of extract- 
ing the metric via spacelike geodesies can be extended to greater than three bulk 
dimensions, the relation between entanglement entropy and minimal surface area 
dHI) no longer involves the geodesic's proper length, and thus this quantity is no 
longer so readily available from the CFT. Computationally, the two methods (as 
presented here) are of comparable efficiency, although both have scope for further 
optimization. 

Significantly, we demonstrated in section [5] how the two methods can be ap- 
plied together to allow the probing of the most general static, spherically symmetric 
asymptotically AdS spacetimes, with metric of the form of (|26p. This is a notable 
extension to the applicability of either method individually, as whilst part of the 
metric information (i.e. the h(r) of (j26p ) could always be extracted using the space- 
like geodesies, they could never give any insight into f{r). The null geodesies, on 
the other hand, can probe both /(r) and h{r) but without yielding enough informa- 
tion to determine either, without some a priori knowledge of a relationship between 
them. It is the separation of the two functions in the spacelike case, however, which 
allows the methods to be combined so straightforwardly, as one firstly determines 
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an estimate for h{r), then takes this as a known function when analyzing the null 
geodesic data. We concluded by considering the toy model scenario of a gas of radi- 
ation (a "star") in AdS^ and demonstrated how the recovery of the metric allowed 
the pertinent information of the star (its mass and density profiles) to be well es- 
timated. Whilst the estimates produced were accurate down to low r (dependent 
on the choice of step size in both the null and spacelike methods), one had to be 
careful when using the derivative of mfit(r) (the interpolating fit to the mass profile) 
to generate Pfit{r). Although the fit produced was still good, it failed at noticeably 
larger r than the fits for any of the other functions, due to inaccuracies introduced 
via mg^.(r). This could be avoided by using the alternative definition of p(r) in terms 
of f{r), (j37p . provided one first extracted the asymptotic fall off of the density as 
r — > oo. 

Finally, we noted in section [6] that this work can be extended further by con- 
sidering spacetimes which are not spherically symmetric, and by investigating the 
higher dimensional cases where the area of the minimal surface relating to the en- 
tanglement entropy is not simply the length of the corresponding static spacelike 
geodesic. Both avenues have the prospect of yielding fruitful results for metric com- 
putation in AdS/CFT. Also, one could approach the problem of metric extraction 
from a different angle, by investigating whether there are alternative methods avail- 
able which do not involve the use of geodesic probes. If so, it would be interesting 
to see whether these lead to more efficient ways of computing the metric functions 
than those described here. 
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Appendix A 

In section Owe outlined the principle behind the iterative technique of metric extrac- 
tion: determining r^m from the gradient dC/d(f) and then calculating an estimate 
for f{r„iin) at each step by splitting up the integral in ([HH]) and approximating each 
piece separately, beginning the whole process at large r, where the metric is approx- 
imately pure AdS and we can take /(r) « -|- 1. Here we go on to give further 
details of how to set this up, and explicitly write down the equations used in the 
approximation^^. 

^•^The procedure used here is only one of a number of possible methods for discretizing the 
integral; for the purposes of illustrating the principle, this method is both brief and accurate to a 
good degree. 
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Having taken the first step which chooses an T^i large enough so the metric is 
approximately pure AdS, and hence f{rn) = + 1, we can continue as follows. For 
a geodesic with slightly lower angular momentum Jn-i (which can be obtained by 
taking a slightly larger region B on the boundary), we can split up the integral over 
r in ([9]) into two pieces: 

== dr + 2 / == dr (45) 

The first integral in the above can be well approximated by taking a next-to- 
lowest order series expansion about the point r = r„_i (= Jn-i), as the distance 
T~n — 1'n-i is small. For the second integral, we can again use our assumption that 
/(r) = + 1 for r > r„, and overall we obtain for the angular separation of the 
endpoints: 



o AT / rn-rn-i 5/(r„_i) 2r„_i/'(r„_i) / r„ - r^-i ^ 



+ arctan 



rn-if{rn-i) 3\/2 \rn-i f {rn^i) , 



— arctan ( " j (46) 

where we have again taken the limit r^ax ^ '''n > fn-i- Alternatively, one could 
perform similar approximations on the equation for the proper length, (jlOp . to 
obtain: 



+2 log (2 r^a.) - 2 log [jrl-rl^^ + (47) 

In the above expressions we have introduced a further unknown, namely the gradient 
of the function /(r) at the point r = r„_i; this can be eliminated by taking the 
simple linear approximation: 

/'(r„-0 « ^^^"^-^^^"-^^ (48) 

which holds provided the radial distance r„ — r^-i is kept smallEll As we men- 
tioned in section 12. H we can calculate the proper length Cn-i from the relevant 
entanglement entropy expression, in general this is given by: 



'^'^The presence of an f'{r) term deserves further comment: one can avoid introducing it by 
using the lowest order expansion, however, this reduces the overall accuracy of the method. The 
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£ = 4Gg)5A (49) 

where Sa corresponds to the entanglement entropy for the subsystem A in the 
deformed spacetime. Taking the entanglement entropy as a known quantity from 
the CFT, along with the angular separation of the endpoints (which is given simply 
from the length of the subsystem in the CFT) one can calculate the corresponding 
minimum radius r„_i from p7|) . and so our only remaining unknown in both (|46|) 
and (I47p is /(r„_i). We can thus numerically solve either for /(r„_i), and determine 
the metric function at this point. Continuing in a similar fashion, by taking geodesies 
with progressively smaller angular momenta and numerically solving at each step, we 
can iteratively extract the complete metric. For general <j)n-i and Cn-i the integrals 
are split up into (i + 1) pieces; two are approximated as in (ji6]) and (jlTl) . with the 
remaining terms evaluated using Simpson's rule (a polynomial fit to the curve). The 
general expression for can then be written as: 

(pn—i ~ -^n—i + Bn—i + Cn—i (^0) 

where 

3/2 



A — 2\/2 /^"'~*+-^ ffi-i ^fijn—i) ~\~ '^'^n-if' {"^n-i) f 1"n—i+l Tn~i \ (51) 
^n-ifiTn-i) 3\/2 V ^n-i f {fn-i) 

Cn-i = arctan I ^ ^ " | — arctan { ) (52) 



n n—i , 



2r„_,^r4 - (r2_^ - 1) 
are the two approximations we had before, and the Bn-i term is given by: 

T3 Sr^ \J'n-2j+2. - rn-2j+l , , \ ^ a I ^ , / \\ 

Bn-i = 2^ 5 {gn-i{rn-2j+3) + 4: gn-i{rn-2j+2) + 9n-i(^n-2i+l)) 

i=l 

(53) 

for i eveio, and by 



(i-l)/2 

E 



D ('"n-2j+2 — '"n~2j) , , \ , a ( \ , / \\ 

Bn-i = 2^ ^ — [gn-i[rn-2j+2) + 4: gn-i{rn-2j+l) + gn-i{rn-2j)) 

(54) 



detrimental effect of the approximation to the gradient on the accuracy of the estimates is not as 
pronounced as in the method of [TD] due to the use of linear step size in J (and hence r), see section 
3.31 Surprisingly, an alternative integral one might consider when setting up the iteration, which 
allows the higher order series expansion to be used without introducing f'{r), leads to an unstable 
method rather than a more accurate one, see Appendix B. 

^^Using this definition requires a value for the r„+i term, which can be obtained in an identical 
way to that used in determining r-„ 
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for i odd, where we have defined the function 

9n-i{r) = — ^_ — = (55) 

V n — i 

for ease of notation. For the proper length we similarly have that: 

^n—i ~ •^n—i ~l~ ^n—i ~l~ Cn—i (56) 

with 

An-i = 2y/2rn-i\ 77 ^ ^ 



(57) 

C„_i = 2 log (2 w) - 2 log (^r2 - r^,. + ^r^ + 1^ (58) 

i/2 . 2-+3-r 2-+l) 
j=l 

(59) 

for i even, and 

(i-l)/2 , 2-+2-r- 2-) 
Bn~i = ^ ""^•^^^g VdM- (Cn-i(r„_2j+2) + 4 Cn-i(?'n-2i+l) + Cn-i{rn-2j)) 

i=i 

(60) 

for i odd, with the function C, defined by 

Cn-^{r) ^ ] , (61) 



r2 . 

n — i 



Thus we can continue the metric extraction down to r = in the non-singular 
case, or down io r = in the black hole case (see section IT3]l . 

Appendix B 

In the method of the previous appendix, the series expansion we used to approximate 
part of the integral in both ([9]) and (llOp introduced an extra term, f'{r), which we 
then chose to linearly approximate. What appears immediately obvious is that one 
could simply combine the two equations and avoid using any approximation to f'{r) 
at all. Considering the two terms An-i and An-i from ([5T]) and ([57j) respectively, 
we see that: 

A - ~ rn^if/'^ 

y fn—ij \fn-i) 
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and so by considering r„_i (pn-i — ^n-i at each step we eliminate the f'{rn-i) term. 
For completeness, we note that this is equivalent to the formulating the integral as 
follows: beginning with expression ()16p and integrating over J gives: 



which can then be integrated by parts: 

£(J) = J(f){J) - j (l)dJ (64) 
and rewritten using the expression for (p from ([9]): 

C{J) = J (t){J) - I I , drdJ (65) 
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We can now reverse the order of integration, and as the function /(r) has no 
dependence on J, integrate over J . For some specific geodesic with proper length 
Cn-i and angular separation (pn-i on the boundary (to continue with the notation 
from earlier) we thus have that: 



Cn~i=rn-i(t>n-i+ j —^==J I ^ dr (66) 

where we have also used that r^m = J and relabeled the minimum radius as r„_j. 
After splitting up the integral as in Appendix A, the lowest order approximation to 
the integral at r„_j is given by (j62p . and one can then seemingly determine /(r„_j), 
the only unknown, for each i from one to n and hence reconstruct the entire metric 
function f{r). 

Applying this in practice, however, one immediately runs into the same stabil- 
ity problems that occur in the naive approach mentioned in section [3l where one 
attempts to recover both r^m and /(r^m) directly from equations ^ and ([TU]) . 
The method appears inherently unstable to errors, and fails to generate any reliable 
estimate for /(r) at any step size. Interestingly, an almost identical formalism can 
be carried out in the method involving null geodesies (see Appendix C), however, 
unlike in the spacelike case, this method is both stable and highly efficient. Fur- 
ther analysis into what causes the stability/instability of the different methods is 
ongoing. 

Finally, to clarify one further point, we note that the original (naive) method of 
section [3] can be stabilised by introducing a particular regularisation of the proper 
length, where one subtracts off the proper length of a corresponding geodesic in 
pure AdS which probes down to the same depth, rmin- Although this appears to 
not introduce any new information, one should remember that we are working from 
the field theory data, and as such, one does not in fact know the proper length 
of this geodesic, but rather the one which has the same angular separation of the 
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endpoints. Thus using this regularisation is actuahy equivalent to determining the 
minimum radius from (jl7p . using this to determine the length of the corresponding 
geodesic in pure AdS, and then treating r^m as an unknown again in (f50l) and (f56]) . 
This excessive over complication considerably reduces the efficiency of the method, 
as the equations are considerably more complicated to solve for (even numerically) 
at later steps. 



Appendix C 

In section [5] we combine the extraction method presented here with that given in 
|lo|^ to allow metric recovery in the most general static, spherically symmetric 
spacetimes. As the methods are applied sequentially, they require very little modifi- 
cation in order to work together, indeed the spacelike method is only affected by the 
change in notation when we introduce h{r). The method involving null geodesies is 
altered slightly more however, and so is presented in full here. This explicit presen- 
tation also serves to highlight the similarities between the two iterative procedures 
for extracting the metric, which is remarkable given the different origins of the field 
theory data. 

As mentioned in the review in section 14. H we have a relationship between the 
gradient of the endpoints of the null geodesies (see figure [TT]) and the ratio of J to 
E, namely dt/d(j) = y, which can be rewritten as: 

dtjy) d(j){y) 
Integrating over y and then by parts gives: 



t{y)=y<P{y)- j (t>dy (68) 
which can be rewritten by substituting in for (j): 



t{y) = y 4>{y) - I I — ^-7==^=^ dr dy 



f'min 



2 / I y 



f{r) r 



Reversing the order of integration (as the function /(r) has no dependence on 
y) allows us to integrate over y: 



t{y)=ym+ r 2^/W)J^.-K^r (70) 



^''There are two methods for extracting the bulk information proposed in [TO]; here we proceed 
to adapt the second, which is noticeably more efficient in generating estimates for /(r). 
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Thus taking the initial conditions to be (^o, to) = (0, 0), we can say that for any 
endpoint {(j)j,tj) on the boundary we have: 



dt 



(</'j>*j) Jr 



dr 



(71) 



where we hcive rencinied Vfyiifi as Vj . After using the spacehke geodesies to determine 
an estimate for h{r), this then finally becomes: 



dt 



(<l>j,tj) 

which, when coupled with the equation for the minimum r, 



(73) 



allows the metric function /(r) to be reconstructed from the plot of the endpoints, by 
applying a similar iterative method to that described in Appendix A: for the general 
term r„_i, one approximates the integral from r„_i to rn-i+i by the parabolic area 
formula; the integral from r„ to r = oo by taking the spacctime to be pure AdS; 
and the remaining i — 1 integrals by the trapezium rule, to obtain: 



where 



and 



dt 

tn.—i . 



i-1 



(.4^71— i 5^n— i) 

4 



'''n—i ~ -^n—i ~l~ Byi—i + Cn—t 



{fn—i+1 fn—i) ViVn—ii ^n—i+1 , 



(74) 
(75) 



Bn-i = XI - rn-j) {viVn-i, rn-j+l) + ??(yn-i, rn-j)) (76) 



Cn 



2 arctan 



—2 yn-i arctan 



yn- 



(1 - yl-i) rl - yl_ 



where we have defined the function r]{yj,rk) as: 



Viyj,rk) = V hfit (rk)] 



fifk) rl 



(77) 



(78) 
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